Crossover from directed percolation to compact directed percolation 
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We study critical spreading in a surface-modified directed percolation model in which the left- and 
right-most sites have different occupation probabilities than in the bulk. As we vary the probability 
for growth at an edge, the critical exponents switch from the compact directed percolation class 
to ordinary directed percolation. We conclude that the nonuniversality observed in models with 
multiple absorbing configurations cannot be explained as a simple surface effect. 
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Recently, considerable effort has been devoted to un- 
derstanding phase transitions in nonequilibrium systems. 
Many studies of models with a continuous transition to a 
unique absorbing state have established that such tran- 
sitions belong generically to the class of directed perco- 
lation (DP) 01, as predicted by Janssen ||^, Grassberger 
and (for multicomponent models) by Grinstein et al. 
1 1 1 . DP has been of great interest since its introduction 
by Broadbent and Hammersley ||^, and is relevant to a 
vast range of models in biology, chemistry and physics 
1^,0. Examples are catalytic reactions on surfaces 
epidemics transport in porous media pC|, chemical 
reactions ||T^,|l^, self-organized criticality |13|JI^ , elec- 
tric current in a diluted diode network , Reggeon field 
theory [p^|-|l8|, and more recently in damage spreading 
, and models of growing surfaces with a roughening 
transition ||20t] . 

Unlike systems with a unique absorbing configuration, 
understanding of models with multiple absorbing config- 
urations pl[-p5[| is far from complete. In such models 
the critical exponents that govern spreading from a seed 
vary continuously with the particle density in the initial 
confi gur ation, and obey a generalized hyperscaling rela- 
tion |24|. The exponents assume the usual DP values 
only for initial configurations having the "natural" parti- 
cle density — that of absorbing configurations generated 
by the system, running at the critical point. Nonuni- 
versality in these models remains a puzzle. In searching 
for an explanation, and remarking that the exponents 
depend on the environment into which the population 
spreads, one is led to investigate whether nonuniversality 
is also produced by modifying the process at the surface 
of the active region. In this work we consider surface- 
modified DP in one spatial dimension, and find it does 
not show continuously-variable exponents, but rather a 
crossover between two different universality classes, com- 
pact and standard DP. Besides providing a simple ex- 
ample connecting the two classes, our study shows that 
the nonuniversality observed in systems with multiple ab- 
sorbing states is not simply a surface effect. 

In bond directed percolation on the square lattice. 



bonds connect each site {x,t) with (x,t — 1) and {x — 
l,t — 1). Each bond is "wet" with probability p; with 
probability 1 — p a, bond is "dry" . Suppose the origin 
is the only wet site in the layer t ~ 0. One constructs 
the cluster Co connected to the origin using the rule that 
(x, t) € Co iff it is connected, by a wet bond, to (x, t—1) or 
(x — l,t — l) G Cq. For p < Pc such clusters are finite with 
probability 1 ; for p > pc there is a nonzero probability of 
a cluster growing indefinitely. The percolation threshold, 
Pc, marks a continuous phase transition. The model is 
readily generalized by introducing conditional probabil- 
ities P(1|1,0) = P(1|0, 1) = pi and P(l|l,l) = P2 for 
site {x,t) to be wet (=1), given the states of {x,t — I) 
and {x - l,t- I). (P(1|0,0) = 0, naturally.) The re- 
sulting Domany-Kinzel model exhibits a line of critical 
points in the DP class j26|. The endpoint of this line, 
(pi — 1/2, p2 = 1), describes a transition outside the DP 
class; it corresponds instead to compact directed perco- 
lation (GDP) |2^]. The essential difference between DP 
and GDP is that in the latter, transitions from wet to 
dry cannot occur within a string of I's; the evolution of 
a string of I's is governed by a pair of random walks at 
its ends. 
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FIG. 1. The centered square lattice and a growing cluster. 

In the present work we study directed percolation on 

the centered square lattice, in which {x,t) is connected 
to (a; — 1,< — 1), {x,t — I), and (a; + 1,< — 1) by bonds 
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that, as before, are wet with probabihty p (see Fig. 1). 
EarUer studies of this model [^,^, yielded a percola- 
tion threshold of = 0.5387 ± 0.0003. We now modify 
the probabilities for introducing wet sites at the surface 
of the active region. Suppose (x^j, t) is the rightmost wet 
site in layer t. We set the probability for {xr^ + !,< + !) to 
be wet at p' rather than p, and similarly for [xl ^ 1, ^ + 1), 
if {xL,t) is the leftmost wet site in layer t. In Fig. 1, for 
example, the probabilities to introduce wet sites at Ai, 
A2 and A3, are p, p and p', respectively. (Similarly, we 
have probabilities p' , p and p for introducing wet sites at 
Bi, B2 and B3 .) 

We used time- dependent simulations to study the criti- 
cal behavior. The method involves starting with a single 
wet site at t = and following the spread of wet sites 
in a large set of independent realizations. We measured 
the survival probability P{t) (the probability that there 
is at least one wet site at time t), the average number 
of wet sites, n{t), and the average mean square distance 
of spreading from the origin, R^{t). At criticality, these 
quantities follow power laws in the long-time limit, e.g., 
P{t) ~ t^^ , and similarly for the other quantities. In 
the subcritical phase {p < Pc{p')) the survival probabil- 
ity decays exponentially, while in the supercritical phase 
[p > Pc{p')) it approaches Poo = \mit^aoP{t) > 0, the 
ultimate survival probability. If we allow the exponents to 
depend on the surface probability, p' , the simplest gen- 
eralization of the usual scaling hypothesis is 



P{t) - t-'5(p')$(Ati/''ii(p')) 

Ti{t) -t''(p')e(Ati/''ii(p')) 
i?2(t) ^t^(p')rj(Ati/''ii(p')) 



(1) 

(2) 
(3) 



where A = \p — Pcip')\ measures the distance to the 
critical point and i^\\{p') is the critical exponent associ- 
ated with the correlation length in the time direction 
(^11 ~ A~'^ii^P '>). Assuming that the scaling functions $, 
and f2 are nonsingular at the critical point, it follows 



TABLE 1. Values of the critical exponents 
and threshold percolation for different values of 
p'. The values for DP were taken from 

and The numbers in parenthe- 

ses represent the uncertainty in the last digit (s). 
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that P{t), n{t), and n{t) follow pure power laws as 
t ^ 00. A log-log plot of P, fi, or i?^ as a function 
of time should yield a straight line at p = Pc{p'), per- 
mitting one to determine Pc{p') rather precisely. The 
exponents 5, 77, and z are given by the asymptotic slopes 
of the corresponding plots. In addition, the exponent 
j3' controlling the approach to the critical point of the 
survival probability, Poo{p) ^ [p — Pc{p')Y |13 
also depend onp'. Our results (see Fig. 2) show that the 
percolation threshold pc increases with decreasing p' for 
p' < p^^ — 0.53875. This decrease in p' means that the 
percolating cone becomes narrow, and a larger value of 
p is needed for the cluster to survive. 

Table 1 show how the exponents vary with p' . The 
mean number of wet sites in surviving trials is N{t) — 
n{t) I P{t), and the fractal dimension d/ of aggregates sur- 
viving to time t is defined through N{t) R'^f , yielding 
df = 2{S + ri)/z. For DP df = 0.752. In the present 
study df ~ 1 for p' < p^^, which means that for this 
range of p' the clusters are no longer fractal but compact. 
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FIG. 2. Phase diagram in p-p' plane. The two regimes 
(CDP and DP) are shown. 
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FIG. 3. The ultimate survival probabihty. Poo as a func- 
tion of p — Pc(p') in a log-log plot. The inset is a plot of P^o 
versus p. These results were obtained for p' — 0.25 and the 
slope of the line is /3 = 1. 
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FIG. 4. The local slopes —5{t) (upper), r){t) (middle), and 
z{t)/2 (bottom), for the case p' = 0.25. Each picture has five 
curves corresponding from bottom to top p = 0.7506, 0.7508, 
0.7510, 0.7514 and 0.7516. 



From Fig. 3 we see that the exponent controling the 
survival probability (/?' = 0.97 ± 0.04) is very differ- 
ent from that of DP (/?' = /3 = 0.2767 ± 0.0004). For 
p' < p^^, we find /3 = 0, so that the generalized hyper- 
scaling relation, (1 -I- (3/ 13')5 + rj = dz/2, becomes 
d + 7] = dz/2, which describes compact growth The 
data are indeed consistent with compact hyperscaling. 
(For the point nearest Pc, Poo is obtained from data for 
t > 15000, averaged over 10^ independent runs.) Figure 

4 shows the local slopes as a function of inverse time for 
p' — 0.25. Similar curves were obtained for the other 
values of p' . (These results reflect averages over 10^ — 

5 X 10^ independent trials, running to a maximum time 
of 2000 to 6000.) We also checked whether the scaling 
relations, Eqs (1) - (3), are obeyed, by varying to opti- 
mize the data collapse. For p' < p^^ , we observe scaling 
for 1 =2, as expected for CDF. 



FIG. 5. Typical evolutions from a single wet site for 
p = p = 0.45, p = p' = 0.5387 and p = 0.9, p' = 0.1 re- 
spectively. Time increases downward. 

Our numerical results indicate that the model exhibits 
a crossover from directed percolation to compact directed 
percolation as the surface growth rate, p', is varied. We 
can understand this by examining a surface-modified con- 
tact process (CP) Q. The CP (in d space dimensions) 
is a sequentially-updated version of DP (in d + 1 dimen- 
sions). In the CP, each lattice site is either vacant or 
occupied by a particle. Particles die at unit rate, and 
give birth at rate A. (In a birth event a new particle ap- 
pears at one of the sites neighboring the parent, if that 
site is vacant.) The CP exhibits a phase transition in the 
DP class, as A is increased through Ac. 

Now consider the one-dimensional CP, with xr and xl 
denoting the positions of the right- and left-most parti- 
cles at any instant, and modify the rules so that the rate 
for a particle to appear at xr + 1, or at xl — 1, is A'/2. 
If we let L =< xr — xl >, the mean being taken over 
all surviving trials at time t, starting from a single parti- 
cle, then dL/dt = A' — 2g, where g is the mean distance 
from XL to the next occupied site (and similarly at the 
right edge). If we set A' = A, there is a critical point at 
Ac, in which case the interior density approaches zero as 
t oo. Since L ^ t^/^ at the critical point, dL/dt — > 
as i ^ cxD as well. (For A < Ac, dL/dt < 0, and vice- 
versa.) Hence, for A' = A = Ac, 5 must tend to Ac/2 
at long times. The ultimate survival probability Poo can 
only be nonzero if dL/dt > at any finite time. On the 
other hand, (7 is a decreasing function of A and A'. It 
follows that if A' < A = Ac, the process must die. We can 
compensate the reduced surface birth rate by increasing 
A, and should again have survival for g(A, A') < A'/2. 
Having forced A > Ac, surviving trials are compact, i.e., 
have a nonzero particle density as i — s- cx). The region be- 
tween XL and xr is not simply a string of occupied sites, 
as it would be in CDP, but the scaling behavior should 
nonetheless be characteristic of compact growth. When 
the surface growth rate exceeds the bulk critical value, 
we expect a continuous (DP-like) transition at A = Ac. 
For A' > Ac, the advance of the edges may be enhanced 
at short times, but the process cannot survive for A < Ac 
because the interior approaches the vacuum. The den- 
sity in surviving clusters is localized at the edges, so the 
number of particles remains 0(1), and the process dies 
with probability 1. To summarize, our argument shows 
why surface-modified DP exhibits compact growth for 
p > Pc{p') when p' < p^^, and why there is a DP-like 
transition at p^^ for p' > p^^ . As in CDP, the transi- 
tion is first order for p' < p^^ . 

Clearly, the surface probability p' does not play the 
same role as the initial density in models with multiple 
absorbing configurations. Conversely, the continuously- 
variable exponents observed in the latter class of models 
cannot be explained as a simple surface effect. Finally, 
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it is interesting to note that while the surface modifica- 
tion explored in the present work represents a nonlocal 
spatial interaction {xr and Xl are defined globally) , the 
dynamics of models with multiple absorbing configura- 
tions involves long memory effects . 
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